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Abstract Young massive stars produce Far-UV photons which dissociate the molecular 
gas on the surfaces of their parent molecular clouds. Of the many dissociation 
products which result from this "back-reaction", atomic hydrogen H I is one of 
the easiest to observe through its radio 21-cm hyperfine line emission. In this 
paper I first review the physics of this process and describe a simplified model 
which has been developed to permit an approximate computation of the col- 
umn density of photodissociated H I which appears on the surfaces of molecular 
clouds. I then review several features of the H I morphology of galaxies on a 
variety of length scales and describe how photodissociation might account for 
some of these observations. Finally, I discuss several consequences which fol- 
low if this view of the origin of HI in galaxies continues to be successful. 
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1. Introduction 

I want to begin this talk with a brief discussion of some aspects of the astro- 
physics of photodissociation regions (PDRs) in order to make it clear that the 
production of H I from H2 in the ISM is quite inevitable, and that the cycling of 
H 1 <-> H2 is likely to be both continuous and ubiquitous in galaxies. I will then 
move on to a brief description of some of the features of the H 1 morphology 
of galaxies which appear to be amenable to an explanation in terms of PDR 
astrophysics. 

The material presented here is a review of results and views published by 
myself and by others in previous journal and conference papers. New in this 
paper are some preliminary results on the effects of radial gradients in the 
metallicity of the ISM on the overall H 1 distribution in galaxies in the context 
of the photodissociation picture described here. 
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2. H i from photodissociation of H 2 in the ISM 

The physics and astronomy of the photodissociation «-> reformation process 
for molecules in the ISM is a subject of active research, and an excellent review 
of the field was published a few years ago by Hollenbach & Tielens(1999). The 
results have been successfully applied to star-forming regions in the Galaxy on 
length scales of typically 0. 1 - 1 pc. The first indication that photodissociation 
may be operating to affect the large-scale morphology (100-1000 pc) of the 
Hi in galaxies was found in M83 by Allen, Atherton & Tilanus(1986), who 
noticed a clear spatial separation between a particularly well-defined dust lane 
several kiloparsec long and the associated ridge of H I and H n in the southern 
spiral M 83. Appropriately for this conference, M 83 is a barred spiral galaxy, 
and the strong stellar density wave which is apparently driven by the bar has 
produced large streaming velocities in the gas across the arms, leading to a 
measureable separation of various phases of the ISM. 

Dissociation of H 2 by far-UV photons 

It requires a photon of energy > 14.67 eV to lift an H2 molecule directly 
from its ground state to the electronic continuum. Such photons will be rare in 
the ISM since they are strongly absorbed in ionizing H 1 to H II. For this reason 
it was initially thought (Spitzer(1948)) that H2 would be very long-lived in the 
ISM. However, closer examination of the electronic and vibrational energy- 
level diagram for H2 revealed another dissociation channel. In this "fluores- 
cence" process, less energetic photons can raise the H2 molecule to an excited 
electronic state. When the molecule decays to the ground electronic state, it 
can end up in a variety of excited vibrational states. However, vibrational lev- 
els above 14 are so energetic that they break the chemical bond holding the 
H2 molecule together, and two Hi atoms are produced. A quantitative de- 
scription of this process was first given by Stecher & Williams(1967). It starts 
when H2 molecules absorb photons primarily at wavelengths of A ps 110.8 
and ps 100.8 nm through transitions to the to the electronically-excited Ly- 
man -» and Werner -» C 1 !!*) bands. In the sub- 
sequent decay to various vibrationally-excited levels of the ground electronic 
state, ~ 10 — 15% of the H2 molecules will dissociate into two H 1 atoms. Con- 
sidering that even higher electronic states exist in H2, it is clear that the FUV 
spectrum over the whole range from 91.2-1 10.8 nm (13.6-1 1.2 eV) contributes 
to the dissociation. In high-UV-flux environments (Go > 10 4 , Shull(1978)), 
photons with wavelengths as long as ps 185 nm (ps 6.6 eV) can continue to 
create Hi by dissociating "pumped" H2 (X 1 S+ ,2 < v < 14) via addition- 
al Lyman- and Werner-band transitions. Verification that this process actually 
occurs in the ISM was obtained when the predicted UV fluorescence spectrum 
was first observed by Witt et al.(1989) in the Galactic nebula IC 63. 
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Using the notation of Sternberg(1988), the rate at which H2 is dissociated 
by this process in the ISM per unit volume can be written as: 

Rdiss = D X x e-^,1000 x /s(iV2 ) x „ 2j (1) 

in units of H2 molecules dissociated cm -3 sec -1 , where 712 = n(H2), and 

D = the unattenuated H2 photodissociation rate in the average ISRF, 
X = the incident UV intensity scaling factor, 
a = the effective grain absorption cross section per H nucleus 
in the FUV continuum, 
Tgr-.iooo = u(Ni + 2N2) is the dust grain opacity at A ~ 100 nm, and 
n = ni + 2n2 the volume density of H nuclei. 

Since D = 5.43 x 10 -11 s -1 (according to the simplified 3-level model, Stern- 
berg^ 98 8)) the time scale for this process on the surface of a typical GMC 
(p,2 ~ 50 molecules cm -3 ) illuminated by the ISRF is r^ ss = (Dri2) -1 ~ 10 
yr! Note that each dissociation produces two Hi atoms on the cloud surface, 
and the appearance of a layer of H 1 when a FUV radiation field is "switched 
on" (e.g. from the birth of a new O-B star) is instantaneous compared to most 
other time scales in the ISM. 



Formation of H 2 on dust grains 

Formation of H2 occurs in the ISM most efficiently on dust grains (cf. e.g. 
Hollenbach & Tielens(1999) and references cited there). The model for the 
formation rate depends on several parameters (some of which are not accurate- 
ly known), as well as on the nature of the dust grains (which may vary from 
place to place in a galaxy). The usual parametrization is: 

Rforrn = 72 X n X U\ (2) 

in units of H2 molecules cm -3 sec -1 , and n = n\ + 2ri2- The rate coefficient 
for unit density, solar metallicity, and 100K kinetic temperature (roughly the 
ISM in the solar neighborhood) is ~ 1 — 3 x 10~ 17 cm 3 sec -1 . This equation 
is strongly dependent upon the dust-to-gas ratio (S = Ay /Nh) and weakly 
dependent upon the gas temperature (T), since 

72 = 3.0 x 10- 18 (S/S )T^ 2 y F (T) cm 3 a -1 , (3) 

where 5q refers to the value of S in the solar neighborhood, and Uf{T) repre- 
sents the efficiency of H2 formation. The product T 1 / 2 j/i?(T) is thought to be 
constant to within a factor of 2 (Hollenbach et al.(1971)). 

The time scale for this process is r/ orm = (2ri72) -1 ~ 5 x 10 8 /ra yr, and 
this will be the rate-determining time scale in an equilibrium situation where 
photodissociation is balanced by reformation on dust grains. 
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Equilibrium 

In recent years, much effort has gone into calculating the level populations 
of the ro- vibrational lines of the ground state of H2, and the intensities of the 
associated quadrupole line emission. These lines can be observed with space 
missions (e.g. SWAS, ISO) in PDRs in the Galaxy and in the nuclear regions 
of other galaxies. The H 1 column density in a PDR is calculated with the same 
physics used to determine the excitation of the H2 near-infrared fluorescence 
lines. While the computations for those lines are rather complicated, the deter- 
mination of H 1 column can be obtained from a simplified version of the model. 
Furthermore the 21 -cm line emission from Hi is almost always optically thin, 
so the observations usually yield the H 1 column directly for comparison with 
the model. The formation and destruction rates described above are set equal, 
and the equation solved, in the present case for the Hi column density. A 
logarithmic form for the analytic solution to this equation was first given by 
Sternberg(1988); see also Appendix A of this paper for a brief derivation. Oth- 
er relevant references are given in Allen et al.(2004). The model is a simple 
semi-infinite slab geometry in statistical equilibrium with FUV radiation inci- 
dent on one side. The solution gives the steady state H 1 column density along 
a line of sight perpendicular to the face of the slab as a function of x, the in- 
cident UV intensity scaling factor, and the total volume density n of H nuclei. 
Sternberg's result is: 



D = the unattenuated H2 photodissociation rate in the average ISRF, 

72 = the H2 formation rate coefficient on grain surfaces, 

a = the effective grain absorption cross section 

per H nucleus in the FUV continuum, 

X = the incident UV intensity scaling factor, 

N(H 1) = the H 1 column density, 

n = the volume density of H nuclei. 

Equation 4 has been developed using a simplified three-level model for the 
excitation of the H2 molecule and is applicable for low-density (n < 10 4 
cm -3 ), cold (T < 500 K), isothermal, and static conditions, and neglects con- 
tributions to N(H 1) from ion chemistry and direct dissociation by cosmic rays. 
The quantity Q here (not to be confused with Go to be defined momentarily) is 
a dimensionless function of the effective grain absorption cross section a, the 
absorption self-shielding function f s {N2), and the column density of molecu- 
lar hydrogen N2 : 



1 DQ 
1) = — x In x + 1 



a |_72 n 



(4) 



f N 2 




Q = / vf, 



Jo 
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The function Q becomes a constant for large values of N2 due to self-shielding 
(Sternberg(1988)). Using the parameter values in this equation adopted by 
Madden et al.(1993), we have: 

iV(Hi) = 5 x 10 20 x ln[l + (90x/n)] 

where n is in cm~ 3 . This is a steady state model, with H2 continually form- 
ing from H 1 on dust grain surfaces, and H 1 continually forming from H2 by 
photodissociation. 

Equation 4 is strongly dependent upon the dust-to-gas ratio (8 = Ay/Nn) 
and weakly dependent upon the gas temperature (T), since 

a = 1.883 x l(T 21 (<V5o) cm" 2 , 
72 = 3.0 x lO-^OV^/VCO cm 3 s-\ 
Q = (a/a ) 1/2 g O , 

where Sq, gq, and Qq refer to values in the solar neighborhood, and uf{T) 
represents the efficiency of H2 formation. The product T 1 / 2 yi?(T) is thought 
to be constant to within a factor of 2 (Hollenbach et al.(1971)). 

Equation 4 also contains a dependence on the level of obscuration in the 
immediate vicinity of the FUV source. While the variable \ represents the in- 
trinsic FUV flux associated with the star-forming region, we generally observe 
an attenuated FUV flux. Assuming any extinction associated with the star- 
forming region is in the form of an overlying screen of optical depth t(FUV), 
X (observed) = \e - T ( FUV ). However, an accurate correction for this effect 
may be difficult, Since the ISM in the immediate vicinity of the FUV source 
will have been disturbed by stellar winds and any prior supernovae. 

Assuming solar neighborhood values of <to = 1.883 x 10~ 21 cm 2 , D = 
5.43 x 10- 11 s _1 , 720 = 3 x 10~ 17 cm 3 s~\ and Q « 5 x 10" 5 (Stern- 
berg(1988)), and neglecting the weak temperature dependence of 72, equation 
4 becomes: 



n \5{ 



(5) 



where x = x(°bserved)e r ^' 7l/ ). The behavior of TV (Hi) as a function of 
X is displayed for 5/S = 0.2, t(FUV) = and n = 30, n = 300, and 
n = 3000 in Figure 1. These values of <5/<5o and t(FUV) are appropriate 
for the outer regions of M101, as discussed in Smith et al.(2000). Values of 
iV(H 1) > 10 22 cm~ 2 are not likely to be observed as the atomic gas probably 
becomes optically thick at this point: 

/oo 
T s T(v)dv 
-00 
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Figure 1. The Relationship between iV(Hl) and \- The observed and predicted values of 
iV(Hl) are shown as a function of \- The modeled behavior of iV(Hl) assumes values of 
8/Sq = 0.2 and t(FUV) = 0, appropriate for the outer regions of M101. The observations are 
clearly consistent with the physics underlying the photodissociation picture. 

-» 1.82 x 10 18 T s tAi> 
« 10 22 cm~ 2 , 

for spin temperatures of T s ps 100K and profile FWHMs of At; « 20 km s _1 
typical of M101 (Braun(1997)), and for optical depths of r ~ 2.5, correspond- 
ing to a ratio between the brightness and kinetic temperatures of Tb/Tk ~ 0.9. 
This value is appropriate for the highest-brightness regions of M101, as indi- 
cated in Figure 8a of Braun(1997). Figure 1 also shows the measurements for 
each of the 35 candidate PDRs in M101 as analysed by Smith et al.(2000). 
The data indicate that the properties of observed regions in M101 are consis- 
tent with photodissociation of an underlying molecular gas of moderate volume 
density. 

Allen et al.(2004) have recently re-examined equation 5 and compared it 
with the full numerical treatment used in the "standard" Ames model sum- 
marized in Kaufman et al.(1999). A conversion of the FUV flux x use d by 
Sternberg(1988) to the quantity Go used by Kaufman et al.(1999) is first re- 
quired; this is because Sternberg(1988) and Kaufman et al.(1999) use different 
normalisations for the FUV flux (see Appendix B in Allen et al.(2004)). When 
distributed sources illuminate an FUV-opaque PDR over 2ir sr, the conversion 
is x = Gq/0.85 (see Footnote 7 in Hollenbach & Tielens(1999)), resulting in 
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90x/n = 106G /ra. With this change, Allen et al.(2004) fitted the analytic 
expression for AT (Hi) above to the model computations in the range in which 
cosmic ray dissociation is not a major contributor, roughly for Go > 1, n > 10 
cm~ 3 . The result is that no consistent improvement is obtained by using any 
value for the coefficient of Go other than the value 106 deduced above, al- 
though a modest improvement is obtained by using a slightly larger value for 
the leading coefficient in the equation, 7.8 x 10 20 , corresponding to a value of 
1.3 x 10 -21 cm 2 for the effective grain absorption cross section. With these 
small adjustments, the final best-fit equation is: 



iV(Hi) = 



7.8 x 10 20 



In 



106G r 



n 



-1/2 



+ 1 



cm 



(6) 



where n = n(HI) + 2n(i?2) and Go and 5o are the (normalized) FUV flux and 
dust/gas ratio in the ISM of the solar neighborhood. In Figure 2 we show values 
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Figure 2. Contours of constant H I column density iV(H i) in units of 10 21 cm~ 2 in the PDR 
as a function of the density n and incident FUV flux Go for the standard model parameters in 
Kaufman et al.(1999) {dotted lines), and for the analytic approximation of eq. 6 {solid lines). 
The labeled contour values are for the numerical model; the contours for the analytic model 
start from 1CT 3 (lower right corner) and increase to 1CT 2 , 1CT 1 , 1, 3, 5, 7, and end at 10 (upper 
left corner). 

from equation 6 plotted as solid lines together with dotted contour lines from 
the "standard" numerical model for 5 = <5o- The agreement is generally good 
over much of the n-Go parameter space of interest here; differences occur 
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mainly in the top left corner of the diagram, and at low values of FUV flux. 
In the top left corner the analytic formula under-predicts the amount of H I 
column density computed from the standard model by about 30% owing to 
H I production by ion chemistry reactions such as + H2 — > H£ + H, 
HCO+ + e" -» CO + H, and PAH' + H + -> PAH + H (where PAH is a 
polycyclic aromatic hydrocarbon), which are important at high Go and low n. 
At values of Go < 1 the contours of AT (Hi) for the numerical model become 
vertical; this is because the standard model includes a low level of cosmic ray 
ionization which contributes a small amount of H 1 by dissociation of H2 even 
for G = 0. 

It should be noted that both the analytic and the numerical forms depend on 
several rather crude parameters used to describe the properties of the gas - dust 
mixture in the ISM (dust cross section, H 1 sticking coefficients, etc.), and that 
not all proponents of PDR models use the same values for these parameters. 
A coordinated effort is presently taking place among the modellers, first to 
see if different PDR codes can produce the same results when using the same 
numerical values for parameters (no, they can differ, and sometimes by a lot!), 
and second, to try to reach some agreement on an acceptable set of values for 
these parameters. 

There are several noteworthy aspects of this equation: 

■ iV(Hl) depends only on the ratio of Gq/u. Low FUV flux, low densi- 
ty environments in galaxies can produce the same column of H 1 found 
in high flux, high density environments. The difference will be in the 
thickness of the H I layer; much thicker layers of H 1 are associated with 
the low flux, low density environments; 

■ at a given n, 7V(H 1) increases first linearly with Go, but then "saturates" 
and increases only logarithmically after that; 

■ at a given Go, N(H 1) decreases logarithmically with increasing n, and; 

■ the H 1 column decreases as the dust/gas ratio increases. 

3. Time scales 

As discussed in the previous section, the relevant time scale for the produc- 
tion of H 1 on the surfaces of molecular clouds is the formation time for H2 on 
dust grains. The full expression including the dust/gas dependence is: 

5x 10 8 
(n 1 + 2n 2 )x(,5/5o) yrS ' 
Table 1 gives typical values for r 2 in several different environments in the ISM, 
and Table 2 lists a number of time scales set by other processes in galaxy disks. 
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Table 1. H I production time scales in typical ISM environments. 



Environment 


n 


S/S 


r (years) 


Solar metallicity GMC 


100 


1 


~ 3 x 10 6 


Intercloud gas 


10 


1 


~ 3 x 10 7 


Outer galaxy GMC 


100 


0.1 


~ 3 x 10 7 



Table 2. Other time scales in galaxies. 



Situation Time scale (years) 

GMC crossing time R/AV ~ 7 x 10 6 

Spiral arm crossing time ~ 5 x 10 7 

B3 star lifetime for FUV production ~5x 10 7 

Galaxy rotation time ~5x 10 8 

Hubble time ~ 10 10 



It is clear that the time scale for H I production by photodissociation of H2 
is short enough to be relevant. For instance, H 1 is produced in the same time 
scales as that of molecular cloud formation (Pringle, Allen & Lubow(2001)) 
and of the star formation in those molecular clouds (Elmegreen(2000)), and in 
a tenth of the lifetime of a typical B star for FUV photon producation. This 
must be the reason why we see H I in regions where young stars form, and 
why H 1 appears in "rims" around clusters of B stars in spiral arms. H I also 
appears in a time short compared to the time it takes for a GMC to cross a spiral 
arm, and the H 1 will correspondingly disappear as the FUV production ceases 
further down stream and the gas reverts to a predominantly molecular form. 
Even in the far outer parts of galaxies we see that the H 1 time scale is still only 
about 10% of the rotation time of the galaxy, so for undisturbed galaxies we 
can expect an equilibrium to obtain between H 1 and H2 even in these sparse 
environments. 

4. Major features in the H i morphology of galaxies 

Our picture of the major features of the H 1 morphology in disk galaxies has 
closely followed the steady improvements in angular resolution of centimeter- 
wave radio telescopes, first with filled apertures (Dwingeloo, NRAO 300', 
Parkes, Arecibo, GBT...), and later with synthesis imaging instruments (West- 
erbork, VLA, ATNF...). 
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Figure 3. Radial distribution of Hi surface brightness for the nearby giant Sc galaxy 
NGC5457 = M 101 obtained by averaging the Hi data in annular elliptical rings. From 
Braun(1997), adjusted to an assumed distance of 5.4 Mpc. R25 for this galaxy is 13.5' w 21 
kpc. 

Radial distribution of H i surface brightness 

A typical averaged H 1 radial surface brightness profile of a nearby giant Sc 
galaxy is shown in Figure 3. The main features to notice are: 

■ the central depression; 

■ the "flat top", typically at a level of 5 — 10 x 10 20 cm~ 2 , and; 

■ the long "tail" to faint levels in the distant outer parts. 

What parts of this plot can be explained by the photodissociation picture 
descibed in the previous section? An explanation for the "flat top" was actually 
suggested in terms of the photodissociation picture nearly 20 years ago by 
Shaya & Federman(1987), but I fear that paper has been widely ignored by 
most workers in the field of extragalactic H 1 and H2. A good look at equation 
6 makes it clear: owing to the logarithmic dependence of N(H 1) on the ratio 
Go/n and the fact that the largest concentrations of young stars also go along 
with regions of highest gas density, we ought not to expect values of H 1 column 
density much in excess of a few times the coefficient in front of the log. So 
even in the most active starbursting regions of galaxies, we are not likely to 
observe values of the Hi column density much in excess of ~ few x 10 21 
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cm -2 . As to the decline in the outer regions, this is likely to be a combined 
effect of a declining ambient FUV flux and a declining area filling factor, since 
Smith et al.(2000) have shown that the total gas density n on the surfaces of 
GMCs located in the neighborhood of massive young stars does not appear to 
change much with radius. 

Finally, the depression in the inner parts can be explained as a consequence 
of the general increase in the metallicity of the ISM in the inner parts of galaxy 
disks. In Figure 4 my student colleague Ben Waghorn has fitted equation 6 to 
the combined radial data on the FUV distribution and the metallicity gradient 
in M 101. The free parameters are the Hi area filling factor (taken here to be 
0.3 everywhere), and the GMC gas volume density n (fits shown for 100, 200, 
and 300 cm -3 ). We see that, in spite of the strong increase in FUV flux Go 
in the inner parts of the galaxy, the rapid rise of the dust/gas ratio (assumed 
proportional to the O/H ratio) actually results in a decrease in N(H i), as ob- 
served, and the quantitative fit is also reasonable. I note here that this result has 
already been described by Smith et al.(2000) for a small subset of young star 
clusters in M 101 accounting for only a few percent of the total Hi content of 
the galaxy; what we are now seeing is that the same explanation appears viable 
for all the H I in the galaxy. We have examined nearly a dozen nearby spirals 
in this way, and find reasonable agreement for about half of them. The other 
half show an indication that there is more H I present than FUV-related pho- 
todissociation can explain. Interestingly, these galaxies nearly all have bright 
nonthermal radio continuum disks, suggesting that there is a component of the 
H I being maintained from dissociation by cosmic rays penetrating throughout 
the GMCs. This work is ongoing. 

Spiral structure 

In the Introduction to this paper I pointed out that it was thanks to a strong 
density wave in the southern barred spiral M 83 that the importance of pho- 
todissociation in affecting the morphology of galaxies on the large scale was 
first unmasked. Other studies have followed on M83 and on other galaxies 
(M 51, M 100; see Smith et al.(2000) for references) and have generally agreed 
that the initial interpretation in terms of photodissociation remains the most 
viable option. The separation in the case of M 83 is about 250 pc, and aris- 
es because of the difference between the spiral pattern speed and the rotation 
speed of the gas, coupled with the time for collapse of GMCs and the time that 
a massive young star lives on the main sequence. 

The clear existence of spiral features in the H I distribution of a galaxy ex- 
ternal to our own was perhaps first convincingly demonstrated for M 101 by 
Allen, Goss, & Van Woerden(1973). Although this galaxy apparently does not 
have a very strong density wave, the Hi is arranged in thin spiral segments 
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Figure 4. Hi column density (top), FUV brightness (middle), and 0/H ratio (lower) as a 
function of radius for M 101. Equation 6 has been used to produce the fits shown in the top 
panel for an area filling factor of 0.3 and n = 100, 200, & 300 cm . From work in progress 
by Allen, Waghom, & Heiner. 

which appear in the inner disk and can be traced over a large part of the main 
body of the galaxy right out to beyond i?25- Figure 5 shows the Hi image 
(grey, kindly provided in digital form by R. Braun) with the FUV contours su- 
perposed (Smith et al.(2000)). The correspondence is excellent, at least as far 
out in the disk as the UIT data extend. We expect to see the FUV image and 
spiral features grow further when the GALEX data become available later this 
year, and we can confidently predict that the close correspondence with the H I 
will continue. 

H i arcs and blisters 

The first study to successfully identify the characteristic PDR "arc" or "blan- 
ket" morphology of H I in close association with far-UV sources in a nearby ga- 
laxy was carried out on M81 by Allen et al.(1997). The problem is, of course, 
to obtain sufficient linear resolution (~ 100 pc) in the H I observations to per- 
mit one to identify the morphology of the PDR structures. An important point 
to note is that the best "correlation" is between the H I and the far-UV, and not 
between the H I and the Ha. 

A study similar to that done on M81 but with more quantitative results has 
been carried out on M101 by Smith et al.(2000), who used VLA-Hi and UIT 



Photodissociation and the Morphology of Hi in Galaxies 



13 



This figure is available as fig5.jpg 



Figure 5. The distribution of atomic hydrogen as seen at a resolution of 6" « 200 pc in 
M 101 (Braun(1995)). Thin, prominent spiral arms wind outwards to beyond the optical disk at 
i?25 ~ 13.5', with the outer Hi arms generally traceable back to arms which start in the main 
body of the disk. Contours of the Far-UV emission recorded by UIT are superposed. From the 
PDR analysis by Smith et al.(2000). 



far-UV data to identify and measure PDRs over the whole extent of the M101 
disk. From these observations they derived the volume density of the H2 in 
the adjacent GMCs in the context of the PDR model. Figure 6 shows the best 
estimate of the H2 volume densities of GMCs near a sample of 35 young star 
clusters. The range in density (30 - 1000 cm" 3 ) is typical for GMCs in our 
Galaxy, lending support to the use of the PDR picture, and also shows little 
trend with galactocentric distance. 

It must be mentioned here that Braun(1997) has offered a different interpre- 
tation of the discrete H i-bright features, which he called the "High-Brightness 
Network" and identified with the "Cold Neutral Medium" phase of the two- 
phase model for the ISM (Field, Goldsmith, & Habing(1969), see also Wolfrre 
et al.(1995) for a more recent discussion). However, in a recent paper, Wolfrre 
et al.(2003) favor the interpretation of Smith et al.(2000) in terms of PDR- 
generated H I . 

There is other, IR spectral evidence that PDRs are important for understand- 
ing the physics of the ISM in galaxy disks. KAO observations of the 158/mi 
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Figure 6. Total gas volume density in GMCs near a sample of 35 young star clusters in M 101. 
This is all H2 deep within the cloud. See Figure 19 in Smith et al.(2000) for further details. 

C 11 line suggest that as much as 70%-80% of the H 1 in NGC 6946 could be 
produced by photodissociation (Madden et al.(1993)), and ISO spectra in the 
mid-IR indicate that the bulk of the mid-IR emission from galaxy disks arises 
in PDRs (Laurent et al.(1999); Vigroux et al.(1999)). 

5. Some implications 

So, you say, let's suppose I am right in my view that Hi is not the fuel for 
the star formation process in galaxies, but merely the smoke from it! Why 
does this matter? Well, the idea that H 1 directly and quantitatively tracks the 
main component of the ISM in galaxy disks is a basic tenet of our current view 
of star formation on the large scale in galaxies, and it may take a few moments 
to consider the alternatives and the consequences of "shifting the paradigm". I 
can think of at least 3 consequences at the moment: 

■ There must be significantly more gas present in galaxies in the form of 
"cold" H2. The H 1 is showing us mainly only the surfaces of molecular 
clouds, and the PDR model by itself does not provide a prescription for 
how to go from the H 1 as a surface phenomena to the H2 in the volume of 
the clouds. But we can confidently predict that more H2 will be present 
than we currently think. H2 in amounts from 2 - 5 times that of the 



-4 



Photodissociation and the Morphology of Hi in Galaxies 



15 



known H I could probably be "hiding" in galaxy disks without clearly 
violating any known constraints. 

■ The far outer parts of galaxy disks, where small amounts of H I appear 
with only sparse star formation, are prime sites for "hiding" H2. A fo- 
cussed effort to find such gas ought to be made. Possibilities for detec- 
tion may include measurements of dust opacity (it is likely to be too cold 
to emit any appreciable amounts of Far-IR continuum or line emission, 
but this needs to be considered carefully 1 ), and molecular absorption 
lines. The anomalous absorption of the Cosmic Microwave Background 
by the 6 and 2-cm lines of formaldehyde H2CO are intriguing possibili- 
ties which ought to be explored further. 

■ We are interpreting the "Schmidt Law" for global star formation back- 
wards. This is a particularly far-reaching consequence of the photodis- 
sociation picture favored here. The situation has been described by 
Allen(2002). Basically, we have inverted "cause" and "effect" for many 
years in this discussion, viewing the H 1 column as the cause, and the star 
formation rate (e.g. quantified by the FUV flux) as the effect. The ob- 
served relationship between these two quantities (roughly a power law 
on a log-log plot) is called the "Schmidt Law for Global Star Forma- 
tion", and has been the basis for many attempts to develop a physical 
theory for large-scale star formation in galaxies involving gravitational 
instability in the disk. Such a theory is still incomplete. On the other 
hand, the PDR picture favored here views the FUV flux as the cause and 
the Hi column as the effect, and provides a simple explanation for the 
observed correlation in terms of physics we already know (see Figures 
7a and 7b, from Allen(2002)). 

6. Summary Remarks 

To summarize the latest views on this topic, first a conclusion which has 
been corroborated by several authors and which by now seems quite solid: 

■ The H 1 spiral arms in the inner parts of "grand design" galaxies consist 
mostly (and perhaps even entirely) of photodissociated H2. 

To this I would add the following points established in papers by myself and 
my co-workers: 

■ As well as O stars, B stars born in spiral arms play a major role in this 
process; 

■ The PDR morphological signature is widespread in HI when enough 
linear resolution (< 100 pc) is available, and; 
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Figure 7a. Data showing a correlation 
between average observed 21 -cm line sur- 
face brightness (converted to "H I mass sur- 
face density") on the X-axis and the Far- 
UV surface brightness (taken as a measure 
of the formation rate of massive stars) on 
the Y-axis, from Deharveng et al.(1994). 
"Schmidt Law" fits are shown, with index- 
es of 1 (dashed line) and 0.6 (solid line). 



Figure 7b. The data plotted with axes in- 
verted so as to emphasize the explanation 
in terms of photodissociation. The solid 
curves are the models of H I production in 
PDRs described briefly in the text, and are 
labelled with the proton volume densities 
of the parent GMCs. The Hi area filling 
factor is assumed to be 0.30 over the disk 
of the galaxy. 



■ since there seems to be no reason to have more than one H I formation 
mechanism, I conclude that both the inner and the (far) outer H I arms in 
spirals are photodissociated H2. 

New results on the radial distributions of H 1 described in this review are 
tantalizing, but need further work: 

■ The general shape of the H 1 radial distributions in many spirals appears 
to be amenable to explanation in the context of a simple photodissocia- 
tion model: 

- variations in H2 density, FUV intensity, and dust/gas ratio can con- 
trol the appearance of H 1, and; 

- Additional H 1 may be produced by an elevated flux of cosmic rays 
in some galaxies. 

Photodissociation of H2 can explain a number of features in the morphology 
of Hi in galaxies on scales from 100 pc to 10 kpc. Note that this process is 
bound to produce H 1 as long as H2 and FUV photons are present; the real 
question we need to answer is "how much?", i.e. what fraction of the total 
H 1 content in a galaxy is cycling repeatedly through a molecular phase on 
time scales which are short compared to the dynamical evolution time of the 
galaxy? If this fraction proves to be large, then there must be an even larger 
reservoir of H2 present to sustain it. Estimating the quantities of cool/cold H2 
hiding in the ISM of disk galaxies is pure guesswork at the moment, but factors 
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of 2 - 5 times that in the form of H I may not be unreasonable, and with a spatial 
distribution as extended as the H I itself. 
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Appendix: Derivation of the H i column density in equation 4 

Equation 4 is obtained by setting the rate of formation of H2 on grains equal to the rate of 
destruction through photodissociation by far-UV photons. Using the notation of §2: 

Rform X Tl X m = Rdiss X 7l 2 (A. 1) 

= J)Goxe^ lm x/ s (JV 2 )xn2. 

For a simple ID slab geometry at constant density n = m + 2n 2 , we can write m = dNi /dx, 
etc., so that the equilibrium equation becomes: 

Rform x n x {dN 1 /dx) = DGo x e ~"( N i+ 2N 2) x j s (jv 2 ) x (dN 2 /dx). (A.2) 

This can be written as a simple first-order separable differential equation in the column densities: 

e +aNl x dJVi = ^/ s (iV2)e- 2CTiV2 x dN 2 . (A.3) 
Rn 

This can be integrated through the entire slab on the half-plane x > to give: 



c \x— oo c \x— 



/ f s (N 2 )e- 2aN2 adN 2 . (A.4) 
Jo 



Rn 

The integral over N 2 on the RHS is just some number, call it G, so that: 

Notes 

1. The paper by Boulanger at this meeting has made a start in this regard. 

References 

Allen, RJ. 2002, in Seeing Through the Dust, eds. A.R. Taylor, TL. Landecker, & A.G. Willis 

(ASP Conference Series, Vol. 276), 288 
Allen, R.J., Goss, W.M, & Van Woerden, H. 1973, A&A, 29, 447 
Allen, R.J., Atherton, P. D., & Tilanus, R. P. J. 1986, Nature, 319, 296 
Allen, R.J., Knapen, J. H., Bohlin, R., & Stecher, T. P. 1997, ApJ, 487, 171 
Allen, R.J., Heaton, H.I., & Kaufman, M.J. 2004, ApJ, 608, 314. 
Braun, R. 1995, A&AS, 114, 409 



18 



Braun, R. 1997, ApJ, 484, 637 

Deharveng, J.-M., Sasseen, T.P., Buat, V., Bowyer, S., Lampton, M., & Wu, X. 1994, A&A, 
289,715 

Elmegreen, B.G. 2000, ApJ, 530, 277 

Field, G.B., Goldsmith, D.W., & Habing, H.J. 1969, ApJ, 155, L149 

Hollenbach, D.J., Werner, M.W., & Salpeter, E.E. 1971, ApJ, 163, 165 

Hollenbach, D.J., & Tielens, A.G.G.M. 1999, Revs. Mod. Phys. 71, 173 

Kaufman, M.J., Wolfire, M.G., Hollenbach, D.J., & Luhman, M.L. 1999, ApJ, 527, 795 

Laurent, O., Mirabel, I. E, Charmandaris, V., Gallais, P., Vigroux, L., & Cesarsky, C. J. 1999, 

in The Universe as seen by ISO, ed.P. Cox & ME. Kessler (Noordwijk; ESA), 913 
Madden, S.C., Geis, N., Genzel, R., Herrmann, E, Jackson, J., Poglitsch, A., Stacey, G.J., & 

Townes, C.H. 1993, ApJ, 407, 579 
Pringle, J.E., Allen, R.J., & Lubow, S.H. 2001, MNRAS, 327, 663 
Shaya, E.J., & Federman, S.R. 1987, ApJ, 319, 76 
Shull, J.M. 1978, ApJ, 219, 877 

Smith, DA., Allen, R.J., Bohlin, R.C., Nicholson, N, & Stecher, T.P. 2000, ApJ, 538, 608 

Spitzer, L., Jr. 1948, ApJ, 107, 6 

Stecher, T.P., & Williams, DA. 1967, ApJ, 149, L29 

Sternberg, A. 1988, ApJ, 332, 400 

Vigroux, L., Charmandaris, P., Gallais, P., Laurent, O., Madden, S., Mirabel, E, Roussel, H, 
Sauvage, M., & Tran, D. 1999, in The Universe as seen by ISO, ed.P. Cox & ME. Kessler 
(Noordwijk; ESA), 805 
Witt, A.N., Stecher, T.P, Boroson, T.A., & Bohlin, R.C. 1989, ApJ, 336, L21 
Wolfire, M.G., Hollenbach, D., McKee, C.E, Tielens, A.G.G.M., & Bakes, E.L.O. 1995, ApJ, 
443, 152 

Wolfire, M.G, McKee, C.E, Hollenbach, D., & Tielens, A.G.G.M. 2003, ApJ, 587, 278 



This figure "fig5.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0407004vl 



